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Abstract 

The procedure to calculate masses and matrix-elements in the presence of mixing 
of the basis states is explained in detail. We then apply this procedure to the two- 
loop calculation in Chiral Perturbation Theory of pseudoscalar masses and decay 
constants including quark mass isospin breaking. These results are used to update 
our analysis of done previously and obtain a value of m u /md in addition to values 
for the low-energy-constants L\. 



PACS numbers: 12.15.Ff,13.20.Eb, 11.30.Rd, 12.39.Fe 



1 Introduction 



Experiments in low-energy QCD are becoming more and more precise. It is therefore 
needed that theoretical calculations are performed to similar accuracy. In the purely 
mesonic sector the low-energy effective theory of QCD is known as an effective theory 
in terms of pseudo-scalar mesons only. It is known as Chiral Perturbation Theory (CHPT) 
and was put on a solid theoretical basis in ffl. In the two-flavour sector calculations to 
two- loop order are now customary. In our earlier work H [| , where also references to other 
calculations can be found, we have performed the main calculations in the three flavour 
sector also to two loops in the isospin symmetric approximation. 

One of the remaining uncertainties was the importance of isospin breaking effects. 
Basically in our previous work and many other theoretical calculations isospin breaking 
was only included by guessing its effect and putting it in the uncertainty of the final result. 
An example where we definitely need to go beyond this is in the use of form-factors Q 
to extract an accurate value of the | V us | CKM matrix-element. Here we take a first step in 
that we evaluate within CHPT the strong isospin breaking in masses and decay constants 
to next-to-next-to-leading order. 

This work serves two purposes: it checks the dependence on isospin breaking of the 
determination of the CHPT low-energy constants (LEC's) at 0(p A ) performed in our earlier 
work |3| and allows us to extract information on the quark mass ratio m^/m^. The latter 
is one of the fundamental input parameters in QCD and thus needs to be determined as 
accurately as possible. 

Another motivation is to see how our previous results 0, |J can change with the new 
experimental data on ||. In particular, they contain one more independent measure- 
ment, the slope of the G form- factor, allowing for a cross-check on the CHPT description 
and we will use it to have a first check on the large N c assumptions made in the theoretical 
analysis. 

The plan of the paper is the following : We describe in Sect. [|the changes needed in the 
presence of mixing of the external states. We first discuss it for the masses and second for 
matrix-elements. We then sketch all the computed quantities, Sect. [|, and follow it with a 
discussion of the electromagnetic corrections in Sect. The latter remain a sizable source 
of uncertainty in our results. In Sect. § we list all the assumptions on which the calculation 
relies. Sect, ^introduces one of the main inputs, the form-factors, discussing the main 
differences in the data and parametrization used in this analysis w.r.t. the previous one. 
Sect. |7| presents a brief summary of our fitting procedure with all the inputs and output 
variables. Sect. ^| is devoted to the update of the results presented already in 0, |3J. In 
Sect. |] we discuss the impact of the outputs on the quark mass ratios. And finally we 
briefly summarize our findings in Sect. |T0| . 

The formulas are of such a length that they cannot be presented in a manuscript of 
reasonable length and we have therefore not included them. For an introduction to CHPT 
we refer to ||. 
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Figure 1: The diagrams contributing to the two-point function Gij. The filled circles 
correspond to one-particle-irreducible diagrams and the lines to lowest order propagators. 

2 Formalism 

In the presence of mixing, where each external field can couple to more than one-particle 
state, we have to generalize the discussion of masses and amplitudes given in || somewhat. 

We first discuss how the masses can be obtained from a general two-point function 
and later the generalization of the wave-function-renormalization part for matrix-element 
calculations. The notation will be chosen appropriate for the tt -t] mixing case, where the 
inverse of a two-by-two matrix can be written explicitly in a simple manner. 

2.1 Masses 

In terms of the lowest order fields <pi we define the two-point function 

Gij(p 2 , m? , F ) = J d d xe ip - x (U\T{Ux)H®)} |0>. (1) 

We have suppressed the dependence on everything except p 2 , F Q and the lowest order 
masses m i0 . The diagrams contributing to this are depicted in Fig. [I]. Writing the lowest 
order propagator as zPy and the one-particle-irreducible contributions as iUki, pictorially 
illki = k ~0~ l , we can obtain the full two-point function in matrix notation as 

G = iP + iP(miP) + iP(miP) 2 + . . . = iP (1 + TIP)' 1 . (2) 

In the case of no mixing, II, P and G are all diagonal matrices and the solution is simple. 
The masses of the particles correspond to the poles which in the diagonal case are given 
by the zeros of p 2 — itlq — II as described in 0. In the general case they are again given 
by the poles in G as a matrix. These can be most easily found as places where the inverse 
of G has a zero eigenvalue or the determinant vanishes 

det(G- 1 ) = det ((1 + IIP) P' 1 ) = det (P" 1 + n) = . (3) 

It is this equation we now like to solve perturbatively in our case of tt -i] mixing. The 
isospin eigenstates used in the Lagrangian are denoted by and r)$. To lowest order, we 
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can diagonalize G and thus by definition P by a simple rotation 

7T° = cos(e) tt° + sin(e) r] s , (4) 
r] = — sin(e) + cos(e) r^g . 

where the lowest order mixing angle e satisfies 

tan2e=^ " ^, (5) 

2 m s — m 

with m = (m u + m d )/2, the average of the lightest quark masses. In the basis defined by 
7T°, r], the Kaon and charged pion fields, we have 

Pa = (6) 

and IT starts by definition only at next-to-leading order, 0(p A ), but does not need to be 
diagonal. The lowest order mixing angle e does not need to be small, we will keep it to all 
orders throughout our calculations. 

With the previous prescriptions P is diagonal and II is blockdiagonal enforced by the 
various other conserved quantities like, for instance, parity. Then each block can be treated 
separately. The charged pion and Kaon do obey a simple diagonal equation without mixing, 
therefore only the 7r°-r) block needs special consideration. 

In order to obtain the chiral expansion of the masses, we need to find the zeros of the 
determinant in Eq. (|3|) function of 

p 2 = ml hys = ml + ( m *r + (m^ + ..., (7) 

where the superindices refer to the chiral order. In the same way we expand the one- 
particle-irreducible self interaction, II (m^ , mf , F ), in chiral orders 

n ( m jL/ S , F o) = n^tm^s, m%, F Q ) + U i6) {m 2 phys , m 2 i0 , F ) + . . . . (8) 

Using 3 and 8 as labels for the 7r° and rj fields defined in Eq. (f|) and temporarily suppressing 
the functional dependence in II the equation becomes to 0(p 6 ) 

{ml hys - ml, + 1$ + 1$) [ml hys - m^ + 1$ + 1$) - (ng + U^f = , (9) 

where the matrix II is symmetric. We insert at this point the expansion of vn? vhys and then 
either the lowest order mass is m| or mg . For definiteness hereafter only the pion mass 
case will be considered. Thus we obtain 

ml = m 2 3Q , 

Kf) = -ngVLm^Fo), 



[m 2 r = -n$( m l , m l,F )-(m^ ^nSV.mS,, 



-m 
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0_ 

dp 2 

U^(ml ,m%,F )Y . (10) 



Notice that we have chosen here to use the strict expansion where we always express the 
arguments of the selfenergy in lowest order quantities. For numerical results later we will 
express all masses in terms of their physical masses and Fq in terms of F n ± and reorder 
the series accordingly. 

We have only used the fact here that a lowest order mixing angle is well-defined. There 
is of course no reason for this to be true at higher orders and indeed it is not the case 
already at order p A [p]]. When photonic corrections are included, the possibility of mixing 
of states with additional soft photons of course introduces the usual infrared problems and 
provides another reason for not having a simple mixing angle for the physical states 0. 

Another alternative is to make a second, p 2 -dependent, field redefinition such that 
P^ 1 + Tl\f (p 2 , m 2 , F ), (i 7^ j) vanishes. In principle in terms of those fields the expressions 
for the matrix-elements derived below are somewhat simpler. The additional complications 
in calculating the diagrams is the reason we work with fields that are only diagonalized to 
lowest order. 




2.2 Matrix-elements 

To calculate any observable we follow the usual LSZ reduction formalism. The matrix- 
element in momentum space for any n incoming or outgoing states is 

A H ... ln = [ 7 =fc£= 1 n - ml) G h .. An (h, ...,k n ). (11) 

i = l i i 

The function (fci, . . . , k n ) is the full n-point Green function generated by the n fields 

fa^ki), . . . , fa n (k n ). The coefficients Zi are defined via 

Gu( P 2 « m 2 phys ) = - Z \ + 0(1) (12) 
p 2 - mf phys 

and are often referred to as wave- function- renormalization. In the diagonal case this is 
the procedure described in detail in ||. In the case where mixing is present, the formulas 
become more cumbersome. We will for simplicity discuss only the case of one external 
leg but the generalization should be obvious. In terms of the Green-function Qj which is 
one-particle-irreducible and amputated on the j-leg which can connect to the fa external 
state we then get 

A = -7w J^M - ™Z)gM)Sj (is) 

with an implicit summation over j. 

The quantity, simply calculable in terms of II, was G^ 1 = —iV, V = P^ 1 + II. Now for 
a two-by-two matrix the inverse is fairly easy 

[P) V 33 (p 2 )Vss(p 2 )-V 2 8 (p 2 ){-V 3S (p 2 ) 7W) J' [ ) 
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Figure 2: A general contribution to G^.,.^. We choose one of the legs, i\ = i, and show 
how the full two-point function Gij contributes on that leg, denoted by the filled circle on 
the right-hand-side. The remainder of the diagrams, the Green-function, reduced on the 
ih leg is depicted by the oval on the left-hand-side and denoted in the text by Qi x i 2 ...i n with 
ih =j- 

where we have indicated the p 2 dependence explicitly. At the physical mass detV(p 2 ) 
vanishes and thus the pole in Eq. (|T2|) appears naturally. For the it and 77 case we obtain 



z ^ ^M^ iKh Z ^ ^M^ H> - (15) 

The final normalized amplitude for the pion case is then 

A 3 = I b -- 1 {^88^3 - P38&} • (16) 

^88^(^33^88-^8) 

all functions are taken at the pion mass. Qj refers to the l.h.s. of Fig. ^|, the dark (red) 
bubble, which contains all the one-particle-irreducible in leg j diagrams. 

The amplitude in Eq. (|T^) contains formally all the chiral orders. Obviously for practical 
purpose one has to restrict to a given order in the chiral series. Therefore we expand V 
and Gi inside Eq. flTE| ) up to 0(p 6 ), but we do not expand these functions in e. 

In our basis we have that (P -1 )^- = Sij(p 2 — m- ) and II only starts at 0(p 4 ) which 
allows to evaluate the expansions explicitly in terms of 

n = nJ 4) + n {6) 



Using the fact that det V = to simplify further the expression, one explicitly obtains for 
the pion case 



J _ C (2) , I C (4) 1 7 (4W2) n 38 g (2) 
A 3 - y 3 + 1 y 3 - 2 Z 33 ^3 - ^2^1 



y 3 ~ 0^33 ^3 



Izgg® + - (zj§) 2 g<?> + z « n « g c2) _ I ( nST\ g m 

2 33 * 3 T g V 33 ) ^ Am2 ^3 g ^ Am2 J ^3 

n ( 4 ) n( 6 ) n^n' 4 ' 1 rK 4 ) 

li -38 MA) 11 38 r (2) il 38 11 88 r {2) J- v (4) ++38 r (2) 



Am 2 8 Am 2 8 Am 2 8 2 66 Am 



where the terms displayed between curly brackets are of 0(p A ) and the ones between 
squared brackets are 0(p 6 ). The formulas are at kf = m%, the physical pion mass and 
Am 2 = (-P _1 ) 88 ( m "iphys) = m K P hys ~ m rjO- We can then rewrite the dependence on the 
lowest order masses and decay constants as before as well in the physical masses and decay 
constants. 

We checked our results in an analytical and numerical way : in terms of the general 
functions sin(e) and cos(e) we find a scale-independent result, while setting e = we recover 
the results of 0. 

3 Analytical results 

In we presented the masses and decay constants and in || in addition the vacuum 
expectation values and the form-factors at next-to-next-to-leading order. Even in the 
latter case we had to resort to numerical approximations for a large part of the expressions 
involving integrals in order to keep formulas to a reasonable size. In the present case, due 
to the presence of isospin breaking effects even the expressions for the masses become very 
cumbersome. In this section we therefore only define the quantities we have calculated and 
how the various chiral orders have been split. 

We have followed as much as possible the philosophy of rewriting the quark masses in 
terms of physical meson masses at 0(p 4 ), with the exception of m u — rrid and the lowest 
order mixing angle e. We have however used Eq. (|5]) as much as possible to reduce the 
number of terms. In addition the lowest order mass relations are extensively used at 0(p 6 ). 

The quantities calculated up to OijP) including isospin breaking due to the quark 
masses to all orders are 

• The 5 meson masses, m^±, rrv^o, ml, m \± an d m 2 K o- 

• The decay constants F n ±, F^± and F^o. 

• The 4-decay constants in the neutral sector. These are the couplings of the physical 
7T° and rj states to the triplet and octet axial current, labeled as F n o 3 , F n o s , F v3 and 
F v8 . In addition we define the quantities e n = F 7r o g /F n o 3 and e v = —F^/F^. 

• The three vacuum expectation values (uu), (dd) and (ss). 



4 Treatment of electromagnetic corrections 

Isospin breaking has two different sources, the quark mass difference m u — irtd and the 
electromagnetic interaction. The former breaking introduces only one new parameter as 
compared to the isospin symmetric analysis of 0, m u /md. The latter can in principle also 
be treated within the framework of CHPT. The corrections to the masses were treated 



in and the framework needed for the decay constants and Ku was worked out in [10|. 



Unfortunately, the general treatment brings in a large number of new unknown constants. 
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Estimating these brings into the game the evaluation of nonleptonic matrix-elements, a 
problem that is beyond our present scope. 

The lowest order is O(e 2 p ) and brings in only one more parameter which can be fixed 
by the difference m^± — m^ ■ The quark mass contribution to this mass difference is very 
small and can be neglected as a good first approximation. This leads to the prediction 



m 2 K ± 



— m 



K° 



cm 



ml 



(19) 



111. At this order there are no corrections 



a result often referred to as Dashen's theorem 
to the decay constants. 

At the next chiral order, 0(e 2 p 2 ), there are corrections to Eq. (|19|). The chiral loga- 
rithms are of the expected size || [12| . The earlier indications of a large contribution also 
from the 0(e 2 p 2 ) LEC's can be found in Jl3|, [TJ]]. These results, but in a less pronounced 
manner, were confirmed by the lattice calculations of |15]. In what follows we will use 



the results of [|T^] which are compatible with the averaged value of [|1^, [H] since it incor- 
porates both a correct long- short- distance-mat chingF] and the chiral logarithms at leading 
and nonleading order in 1/N C . It also agrees with the results of fllS] . The final result of 
16fl was 

(m 2 K± - m 2 K0 ) em = (1.84 ± 0.25) (rnfe ~ m*o) . (20) 
We will use the central value of Eq. (EOf) in our numerical results and Eq. ([19]) to show 



the dependence of the result on this quantity. Numerically the above two correspond to 



(m K± - m K o) em w 2.3 MeV [Eq. (20)] and 1.3 MeV [Eq. (@)] . 



(21) 



We calculate the masses and decay constants with the strong meson masses, i.e. with the 
electromagnetic effects above removed. 

The electromagnetic corrections to the decay constants are only partially treated within 
CHPT. We will use the values as calculated by |18|, |19| as F n ± = 92.4 MeV and F^±jF^± 



where the electromagnetic radiative corrections have been removed. 



The experimental results of [20] were already corrected for radiative corrections. In the 
absence of better calculations and measurements we therefore use our earlier procedure |§ 
again. We calculate this process with the physical mass of the charged Kaon and pion. 



5 Other assumptions 

The calculations in CHPT we have performed are all to 0(p 6 ). We have not used any 
improvement procedure to guess the size of the higher order terms nor can we determine 
experimentally the 0(p 6 ) tree terms. At present the data are simply not constraining 
enough to allow us to do so. The total number of parameters is very large. Thus these 
have to be estimated using a theoretical procedure as discussed in Sect. fO]. In addition, 
the available data and 0(p 6 ) calculations are not sufficient to determine all other inputs, 



the assumptions we use here are given in Sect. p72 . 



1 It is precisely this matching that invalidates the results of p 
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5.1 At next-to-next-to-leading order, 0(p 6 ) 

When the full fits of the 0(p 4 ) LEC's became available |l|] an understanding of their 
general pattern and sizes within about 30% could be obtained using resonance saturation 



2jJ . These results were later refined in the framework of quark-models and large N c [^2 
We will use the former procedure, a few tests at 0(p 6 ) exist, in two-point-functions some 



parameters could be determined [^J and in the pion form-factors two additional tests were 



possible p4 |. All of these were successful within about 30%. The precise estimates we use 



are described in [|. We test this assumption by varying the input in several ways: 

1. Keeping only the two main contributions: vectors and rf. 

2. Varying the input parameters of the resonances by a factor of two. 

3. Choosing different renormalization points where we do the estimates. 

The combination of the above 3 variations results in a rather large variation of the 0(p 6 ) 
constants used. We thus expect that this part does not affect our result within the quoted 
errors. 

5.2 At next-to-leading order, 0(p 4 ) 

Furthermore, using the data presented in p^] we need to fix four more inputs, we choose the 
quark mass ratio m s /rh, and three of the LEC's: L\, L\ and Lg. The new free parameter 
compared to is m M /m^ but now there are two Kaon masses instead of one as fitting 
parameters so m M /m^ can be determined as a function of the inputs. 

For the quark mass ratio m s /m we use as standard value 24. This is a reasonable 



average of the ratio of the sum- rule calculations of rh |^5[ and m s |26[ . We will also check 
the changes of the fits with a rather broad range of inputs for m s /rh G [20, 30]. 

The constant Lg driving the pion charged radius is saturated up to some good extent by 
a Vector Dominance expression. The pole due to the p exchange represents the dominant 
low-energy singularity while the inr and KK cuts generate a tiny correction. This LEC is 
fixed at its standard value of 6.9 x 10~ 3 . Even putting it to zero hardly changes the results 
so this assumption does not bring extra uncertainty. 

The largest uncertainty w.r.t. the masses and decay constants is in the values of L\ 
and Lg. In the large N c limit they vanish and varying the scale fi from 500 MeV to 1 GeV 
changes L\ by about ±0.5 x 10~ 3 and Lg by about ±0.3 x 10~ 3 . A range within zero of 
that size is thus an acceptable range. There are some indications that the values are at 
the end of this region and even somewhat outside for Lg ||27| , |28|. We will make use of two 
experimental data analysis : in the first one we will use (otherwise using p9f barely 
will change the conclusions) and will keep both of them, L\ and Lg at zero for our main fit 
but discussing the variation of the results below. In second place we will use the analysis 
of H , this will allow to remove one of the assumptions for the large N c behaviour of L\ or 
Lg and compare the result with the previous one. Let us note here that there are possibly 
large corrections to the large N c limit in the scalar sector but they are to a large extent 
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due to the chiral loop corrections and are thus included in the present results. Examples 



are n-n scattering || and the pion scalar form-factor ||24j| . A large N c saturation of the 
constants in the latter case provided excellent agreement with experiment. This question 
does however deserve further study. But if we assume that the SU (2) quantities should not 
be too sensitive to m s then we need to assume that L\ and L r G are small. We also obtained 
that the 1/N C suppressed combination (2L\ — U^)jL\ was small providing one more 
test of the large N c assumption. Let us comment briefly about the resonance saturation 
pattern of the LEC. The symmetry breaking part of the effective Lagrangian at 0(p 4 ) 
is determined by L\, . . . , Lg. Something in common to all of these constants is the fact 
that in the low-energy regime none is affected by the exchange of spin 1 resonances but 
of scalars and pseudoscalar particles, sector in which the experimental data are still fuzzy. 
Further more, for L r 4 , ... , Lg the natural scale (Ms — ~ 1 GeV) differs from the one 
of L\,L\ and Lg, set in the last ones by Vector exchange (My ~ M p ~ 0.8 GeV). This 
obviously is reflected in the expansion in powers of m s , some quantities are controlled by 
the parameter (M K /M S ) 2 ~ 1/9 while others present a worse behaviour (M K /M p ) 2 ~ 1/4. 

6 Ku decay parametrization and data 

The experiment E865 was not primarily designed to collect data but it recorded it as 
a background anyway ||. These data were previously analyzed in a way similar to the 



older experiment jZIJ an d presented in |29] . The second and latest analysis || uses a more 
suitable parametrization of the form-factors to fit the data at low-energy |3(| 

F = {f s (s n ) + f eSi }e iS °oM + p- W ave, (22) 
G = {^p(s^) + g e s e }e iS * M + d-wave , 

where K — > n + (p+) n~(pJ) £(pi) v(p v ) decay variables are defined as 

si = (pe + Pv) 2 , s,r = (p+ + P-f ■ (23) 

The analysis of assumes only dependence in the lowest partial wave decomposition, i.e., 
s-wave in F and p-wave in G. Furthermore any sensitivity on si is discarded. As previous 



work showed, 0] , those assumptions where valid in the analysis of |20[ . But in view of the 
errors quoted in || for the form-factor F at threshold, / s (0), this assumption is borderline 
for the new data. 

The form-factors in || where fitted to the following expressions 

F = (fs + f' s s n + f'Jsl)e<^, (24) 
G = (g P + g' p s„)e iSl ^, 

where f s ,f' s , f'J, g p and g' do not depend on any kinematical variable leading to the result 

/ s (0) = 5.75 ±0.02 ±0.08, f' s (0) = 1.06 ± 0.10 ± 0.40 , (25) 
P (O) = 4.66 ±0.47 ±0.07, g'JO) = 0.67 ± 0.10 ± 0.04 . 
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The quadratic slope is rather larger than we expect from our CHPT results 0. We therefore 
use the results of || with a linear fit only 

f s (0) = 5.770 ±0.097, f' a (0) = 0.47 ± 0.15 , (26) 
g p (0) = 4.684 ±0.092, ^(0) = 0.54 ± 0.20 . 

and regards these fits as preliminary till the data are finalized. 



7 The fitting procedure and fit results 



The fitting procedure is the same as the one used in || but slightly adapted to the present 
case. We refer to and the previous sections for a description of the details about extra 
assumptions, besides CHPT, and how electromagnetic corrections were taken into account 
in all the computed physical quantities. Here we only sketch briefly the inputs and outputs 
of the calculations. 

In addition to the lowest order Lagrangian that contains: Bq m u , Bo , Bo m s and 
Fq as free constants (see below) the 0(p A ) Lagrangian is parametrized in terms of 10 un- 
known LEC's not restricted by symmetry: L[, ... , L r lQ . All in all this amounts to 14 free 
parameters in total Q. The constant L r 1Q does not appear to 0(p 6 ) in any of our quan- 
tities, thus remain 13 free parameters to be determined at 0(p % ). Taking into account 
the restrictions of Sect. ||the number of parameters is reduced to 10 or 9 depending if we 
use the results of [|20| or || respectively. Hereafter we will quote the quantities referring 
between parenthesis. The experimental inputs used are four of the five pseudoscalar 
masses, i.e. we discard the use of m^± because its difference with m 2 is purely of electro- 
magnetic origin. Two decay constants F n ± , F K ±. And three (four) inputs from Ku decays, 
Sect. H two form-factors at threshold, / s (0) , g p (0), and one (two) slope f' s (0) , (g' p (0)). In 
all there are nine (ten) inputs and therefore this allows to determine nine (ten) parameters. 
Schematically 



Inputs 



Outputs 



F«±,F K ± 

M 2 ,M 2 K± ,M 2 R0 ,M 2 
F 

Bo rh 
m u /m d 

L\ , L 2 , L 3 , (Li) , L 5 , L 7 , L g 



(27) 



(2* 



Notice that using we obtain between our outputs the order parameter L\. Unfortu- 
nately the Ki4 quantities are not very sensitive to L\ |H] . 



2 Notice that we have one more parameter than in j3[ due to the isospin breaking. 
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Table 1: Results for £[(/f), Fq, Bom and m u /md for the various fits described in the main 
text using [p0 |. Errors are fitting errors as quoted by MINUIT. All L^(fi) values quoted 
have been brought to the scale fi = 0.77 GeV. The first line with m M /m^ (full) is with 
Eq. (pop , while the second line, m u /rrid (simple) is with Dashen's theorem, Eq. fllTf). The 
standard values are m s /m = 24, = 0.336 GeV, S£ = 0, L\ = Lg = and L r 9 = 6.9 10~ 3 . 
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0.0136 


0.0138 


0.0158 


0.0110 


0.0136 


(5m,, [MeV] 


0.32 


0.27 


0.16 


0.32 


0.32 


0.23 


-0.07 


0.63 


0.36 


m u /m d (full) 


0.46±0.05 


0.52 


0.43 


0.45 


0.46 


0.46 


0.47 


0.33 


0.46 


m u /md (simple) 


0.52±0.05 


0.58 


0.50 


0.53 


0.53 


0.53 


0.54 


0.43 


0.53 


changed 




o( P A ) 


m s /m 


Vs~e 




T r . j r 

^4; ^6 


1* 


A* 


5(0) 


quantity 






26 


0.1 


0.293 


-0.3; -0.2 


0.5 


1.0 


4.93 


Unit 








GeV 


GeV 


10~ 3 


GeV 


GeV 





The errors used in the fit were the experimental errors on f s , f' s and g p , (g' p ), 0.01 for 
Fk/F^ and we required the quark mass ratios determined from the lowest order meson 
masses to be satisfied within 10% accuracy. 

The procedure for the latter was to take the physical meson mass, calculate the 0(p A ) 
and 0(p 6 ) corrections to obtain the lowest order masses of n ± , ir°, iV ± , K° and rj. Since the 
lowest order masses actually were needed in the calculations, they surface in the formulas 
via e and Am 2 , we iterated the procedure till it converged. 

The change in the L£ using |2D] are compared to those reported in for all the 



fits was significantly smaller than the errors quoted there, indicating that the neglect of 
isospin violation was a good approximation^. The changes are always below 0.02 x 10~ 3 
in magnitude. We therefore still consider the main fit of as the standard values for the 
LI 

The remaining electromagnetic effects are estimated by also performing a fit, this time 
with the masses of all the propagating particles inside the loop set to their physical mass 
rather than to the strong mass only. Again, the changes in the results were rather small 
compared to the quoted errors. 

The fit results using [|2D| are presented in Table [TJ This can be compared directly with 
Table 2 in With respect to this we have discarded the comparison with previous works 

3 Note that the change in value is due to an error in the numerical program used in ]3| as described in 
the erratum. 

4 See footnote 0. 
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[!2[] (second column in Table 2 in 0) and fit 6, performed only with vectors and r/ in 
the G\ saturation, otherwise we keep the same notation and only include m u /m d , F and 
Bqitl as output parameters. We also performed the fit with varying the resonance input 
parameters by a factor of two with similar results. 

Let us summarize our findings in Table [I]: In the Main fit we use a strict standard 
approach, considering the central value of the Zweig Rule violating term, L\ ~ Lg = 0, 
r = 24 and the scale is set to M p . The second column compares each quantity with the 
corresponding 0(p A ) results using the same input assumptions. Fits 3 and 4 correspond to 
different choices in the kinematical points for the Ki& variables. In fit 5 we allow for a small 
violation of the Zweig Rule. Fits 7 and 8 show different choices of the scale dependence 
ji. Our results are fully /i-independent, but numerically the change reflects our lack of 
knowledge about the scale where the resonance saturation for the LEC's, C^, is valid. 



Finally fit 9 presents a weighted average of the p-wave form-factor g p (0) |2^. |33| . 

There are no changes w.r.t. the numbers in HQ. Our results are therefore stable against 
the isospin breaking effects. Furthermore the values of L\ and Lg are compatible with the 
recent lattice simulations [33 ]. 



In order to show the changes introduced by the analysis of we also include Table [|. 
Notice that the values quoted there are only indicative and should be taken as preliminary 
till the final E865 analysis is available. As explained in Section |6| we use the linear fit of 
H to perform the fit. 

In Table ^] we have displayed the following results : fit 10 is the direct equivalent of the 
main fit, and the next column are the corresponding 0(p 4 ) values. Fit 11 is the equivalent 
of the main fit but leaving L\ as a free parameter. As one can see there are no significant 
changes w.r.t. assuming L\ ~ and we only obtain a very weak limit on L\. 

This is not conclusive, i.e., one can not claim that Lg = =^> L\ ~ is the only 
possible solution because both constants can be highly correlated. For this reason we show 
in fits 12 and 13 the shift due to different choice of Lg. As one can see from the last three 
columns of Table |2] the values of Lg and L\ follow almost a linear relation. Therefore a big 
deviation from the Zweig Rule value of Lg will probably signal also a deviation for L\. The 
possibility of a small value of Lg and a large value of L\ fixed via the scalar pion radius is 
thus unlikely. 



8 Survey of Applications 

In this section we collect various results in the light of the new values of the main fit in 
Table [I]. We want to stress that we do not study the convergence of CHPT here, we only 
want to check that the total size of corrections is not very large as compared to the lowest 
order. The masses typically have small 0(p 4 ) corrections but sizable 0(p G ) ones. As we 
mentioned in [Q, [| there can be several reasons for this kind of behaviour, we rephrase 
them here for sake of completeness : there is a strong suppression between the LEC's and 

5 See footnote ||. 
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Table 2: Results for Ll(/i), F , B m and m u /md for the various fits described in the main 
text using the inputs of ||. L\ and Lg are given, = means it is an input value. All other 
inputs are as in Table [I]. 





ncru 


U{p ) 


mi i 


nciz 


ni io 




n /i q _i_r> 1 o 


U.oo 


n /i q _i_r> 1 o 


n /i q 


n /i q 
U.4o 


in 3 T r 


n 7Q-un 1 9 


1 

1 .ov 


n 774-n 9D 


u.oo 


D 7D 


10 3 L r 3 


-2.35±0.37 


-2.91 


-2.36±0.40 


-2.40 


-2.32 


10 3 L\ 


=0 


=0 


-0.18±0.85 


-0.46 


0.10 


10 3 Ll 


0.97±0.11 


1.46 


1.08±0.61 


1.07 


1.11 


10 3 U Q 


=0 


=0 


=0 


= -0.2 


= 0.2 


10 3 L T 7 


-0.31±0.14 


-0.49 


-0.35±0.30 


-0.35 


-0.36 


10 3 L£ 


0.60±0.18 


1.00 


0.70±0.59 


0.70 


0.70 


F [MeV] 


87.7 


81.1 


93.9 


105.7 


84.1 


B m [GeV 2 ] 


0.0135 


0.0181 


0.0104 


0.0106 


0.0104 


5m,, [MeV] 


0.28 


0.27 


0.83 


0.53 


1.12 


m u /m d (full) 


0.45±0.05 


0.52 


0.37±0.12 


0.38 


0.38 


iriu/n^d (simple) 


0.52±0.05 


0.58 


0.45±0.27 


0.45 


0.47 



chiral logarithms at 0(p 4 ) and the assumption of resonance saturation using our naive 
scalar picture at 0{p G ) might also be incomplete. 

8.1 Meson masses 

The numerical results for the masses can be presented in various ways. Using the definition 
Eq. (0) we obtain 

my(m^) QCD = 0.742 + 0.007 + 0.251, 
m 2 K ±/{m 2 K ±) QW = 0.693 + 0.027 + 0.280, 

m 2 v /(m 2 v ) phyB = 0.741 - 0.029 + 0.289. (29) 

These results are the same within errors as those given in 0. Hereafter the first quoted 
number refers to the Born approximation, the second to the next-to-leading correction and 
the third to the next-to-next-to-leading contribution. Notice, that as previous results con- 
cerning the masses, even if the 0(p 6 ) contributions are bigger than the 0{p A ) all three ap- 
pear to be of the same size, (m 2 ±/(m 2 ±) QCD ) (6) ~ (rri 2 K± /(m 2 K± ) QCD ) {6) ~ (m 2 / (m^) pliys ) . 
Furthermore, almost all the quantities (besides the 0(p 6 ) LEC's) contributing to the masses 
have the same sign. Thus in order to make 0(p A ) > 0(p 6 ) there must be some fine tuning 
of the 0(p 6 ) LEC's. Our simple not well controlled estimate of the LEC's at 0(p 6 ) can 
thus be the reason for large 0(p 6 ) corrections. 
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The isospin breaking quantities have as expansion 

(m^ -m^±) QCD = (5.09- 0.34 + 1.53) x 1(T 3 GeV 2 , 
(mj ± -m^ cd = (3.96 + 1.54 + 3.17) x 10" 5 GeV 2 . (30) 



A GA ./o = m 2 ri + mlo - - (m) <± + m 2 K o ) - -m 2 ± = . (31) 



QCD 

The Gell-Mann-Okubo relation in the presence of isospin breaking is 

2 / 2 2 \ 2 

3 (™ K ±+m K0 ) - -„. ff 

The lowest order masses of course obey this identically and we use it extensively throughout 
the calculation. The deviation is given by 

A GMO = (0 - 0.0172 - 0.0033) GeV 2 = -0.0205 GeV 2 . (32) 

8.2 Decay constants 

The decay constant have as expansion 

F w ±/F = 1 + 0.135 - 0.075, 
F K ±/F^± = 1 + 0.162 + 0.058, 

i^g/i^i = 1.000 + 0.242 + 0.066 = 1.308, (33) 



and for the isospin violating quantities 



.FLoq — F w ± , . A 
= (1.0 + 0.3 + 0.8) x 10~ 4 



F v ± 



_F^ 
F v s 



0.0141 + 0.0039 - 0.0021 = 0.0159 , 

0.0141 - 0.0010 - 0.0026 = 0.0105 . (34) 



8.3 Vacuum expectation values 

The effective Lagrangian at lowest order contains two order parameters Fq and Bq both 
of them sensitive to the infrared end of the Dirac spectrum [35] but their role in the 
spontaneous chiral symmetry breaking pattern is quite different. The only necessary and 
sufficient condition of SB^S is a nonzero value of the correlator function 

lim % [ d 4 x( n | : UJx)Ri(V>) : | ft) = -~»^> 2 + 0{m q ) , (35) 

where = \{V^ — A^) and = h(V^ + A^) are Noether currents. There is no such 
constraint on B whose size is determined by the vacuum expectation value 

(0\q i q j \0) = -F*B 6 ij O(m q ), (36) 
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and in principle has no lower bound other than vacuum stability, Bq > 0. Even though 
there is no proof available is rather likely that QCD is not realized in the phase B = 



The chiral expansion of the vacuum expectation value ( p6|) can be ordered as 

( | gV | ) = -F%B {1 + ( | q l q j | ) (4) + ( | (fq j | ) (6) + . . .} , (37) 

where the superscripts refer to the chiral order. Although the result is scale independent 
order by order in the CHPT expansion it depends on the QCD renormalization scale Hqcd- 
It follows from the fact that the quark scalar current needs to be defined precisely in QCD. 
Already at 0{p 2 ) B depends on Hqcd- At 0(p 4 ) there appears an additional ambiguity 
via the high energy constant H 2 which is forced by the ambiguity in the subtraction of the 
scalar two-point function 



<0|w|0> = (0|m|0> cMral Hmit - i J d*x(0\ : (qMq)(x)(qq)(0) : |0> + 0(M 2 ) , (38) 

at zero distance. At 0(p 4 ) one can eliminate the constant H 2 and thus obtain a sum rule 
relating the ratios of the isospin asymmetry (0|c&i|0)/(0|mt|0) to the SU(3) asymmetry 
(0|ss|0)/(0|mt|0). In order to quote numbers we fix H 2 assuming scalar dominance 

H' 2 =2L\. (39) 

With this approach we obtain the following expansion 

(0|im|0) = -B Q Fj}[l + 0.271 + 0.106] , 
(0|dd|0) = -B F$[1 + 0.284 + 0.110] , 

(0\ss\0) = -B F 2 [l + 0.964 + 0.289] , (40) 

where the quoted numbers correspond to the lowest order, 0{p 4 ) and 0(p 6 ) respectively. 
The dependence on H 2 can be judged fromQ 

HI = 0, 
(0|««|0> = -B F 2 [1 + 0.264 + 0.104] , 
(0|dd|0) = -B F 2 [1 + 0.271 + 0.107] , 

(0|ss|0) = -£ F 2 [1 + 0.693 + 0.222], (41) 
The vacuum asymmetries are defined as 

(0\3d\0) 



{ | uu | } 



(o 


ss 


0) 


(0 


uu 


0) 



1 - ; '_ ' =l-e s . (42) 



3 The 0(p 6 ) dependence on 7J£ is from rewriting F 2 in in (0|qq|0). 
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Using Partially Conserved Axial Currents (PCAC) these ratios are fixed to be 1, then 
and e s measure the SU (2) and S77(3) breaking of the non-perturbative vacuum respectively. 
We get for them 



e d = -0.0129 - 0.0011 
e s = -0.693 + 0.004, 



(43) 



where the first number corresponds to 0(p 4 ) and the second to 0(p 6 ). As becomes clear 
the breaking inside the 577(2) doublet driven by is very small while the breaking due to 
the presence of m s , e s , is rather huge. There is a plethora of phenomeno logical calculations 
determining Eq. (f|3|) from sum-rules. Their results are very spread and no clear conclusion 
can be reached. Otherwise the results in Eq. (|43|) are in the ball park of those quoted in 
[37] for e d and |38[] for e s while the remaining results |39| are quite off of Eq. (|43"|). 



9 Quark Mass Ratios 

A review of the early quark mass determinations can be found in f40|| . A more recent list 



of references is the quark mass section in the Review of Particle Properties ||T] . 

While sum-rule or lattice techniques allow to extract the value of quark masses via 
some hadronic observable, the use of an effective Lagrangian, as CHPT, only allows to 
determine the relative size of the quark masses, m u ,m d and m s . This approach assumes 
that the quark masses can be treated as a perturbation on the massless QCD Hamiltonian, 
Ho 

7~Lqcd = 7~Lq — J d 3 x [m u uu + m d dd + m s ss] . (44) 

Quark masses depend on the renormalization scale, but if we use a mass-independent 
renormalization scheme like MS they are renormalized multiplicatively and their ratios 
are renormalization-scale-independent. Due to the locality of the symmetries in the chiral 
Lagrangian the quark masses in CHPT can be identified with those in QCD via the Ward 
identities, thus by construction they are the same as those in QCD. 

The values of the relevant LEC's can depend on the precise definition of the scalar 
currents in QCD. In practice, since we have to determine the LEC's from the same data as 



we use to obtain the quark masses some ambiguities might occur |42j| . These ambiguities 



are discussed in where we refer for a more detailed discussion. In the present context 



they are fixed by our large iV c assumptions on the values of L r A and L r 6 . 

To lowest order in the chiral expansion and neglecting electromagnetic corrections we 
have 

ml, -5 = ml ± = B (m u + m d ) , 
m 2 K ± = B (m u + m s ), 
<o = B (m d + m s ), 

= -B (m u + m d + Am s ) + 5 , (45) 
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where 5 is the effect of 7r -r] mixing. The constant of proportionality is determined by the 
quark condensate B = | ( \uu \ ) \/F 2 + 0{m q ). We can use Eq. (|45D to get 



m s _ m 2 K ± + m 2 K0 - m\± 



rh m 2 ± 

TV 

m u m 2 K± - m 2 K0 + m\ ± 
m d m 2 K0 - m 2 K± + m\ ± 



24.2 (real); 25.8 [Eq. (19)]; 25.7 [Eq. ©] , 

0.66 (real); 0.56 [Eq. (||)]; 0.49 [Eq. ©] . (46) 



The three quoted numbers correspond to using the physical masses, electromagnetic cor- 
rections using Dashen's theorem, Eq. (|H]), and e 2 m s corrections included as well, Eq. (|2"0|). 

To 0(p 4 ) we do not have quite enough information anymore to determine both quark 
mass ratios and all the relevant LEC's. To first order in isospin breaking it was noted in 
|1|, that the particular combination 



^ 2 _ m 2 s -fh 2 mj\ 1 + (m u /m d ) = m\ m 2 K - m 2 n 

m 2 d -ml A\ rh 2 J 1 - (m u /m d ) m 2 m 2 K0 - m 2 K± 

could be related to the meson masses without corrections at this chiral order, thus providing 
a tight constraint on the quark mass ratio. As was pointed out in [J2] if the terms rh 2 /m 2 



are discarded Eq. fl47|) leads to an ellipsis where the variables are the ratios m u /md and 
m s /m d . The major semi-axis is given by Q while the minor one is equal to 1. Plugging 
the masses in Eq. fl47j) gives 



Q = 22.0 (21.3) [Eq. (20)]; 24.1 (23.4) [Eq. ([TJD] . (48) 

The numbers are by using the strong meson masses in the last part of Eq. ($7\). The 
presence of higher orders corrections can be seen from the number in brackets, these are 
calculated using the quark mass ratios from the fit with m s /m = 24 at 0(p 4 ) of Table [l] and 
the definition of Q in terms of the quark masses. Again the two sets of numbers correspond 
to using the electromagnetic corrections estimated with e 2 m s corrections included, Eq. (|20|), 
or G(e 2 ) only using Dashen's theorem, Eq. (|T9"D. 

The quantity Q can be determined to the same order from 77 — ^ 37r decays giving 

Q = 22.7 ±0.8. (49) 



The decay rj ^ 3tt is known in CHPT to 0(p ) MM and higher order corrections were 



estimated using dispersion relations [[4I|. In view of the discrepancy between dispersive 
estimates f32| and the full calculation to OijP) fU in the Ku process the error in Eq. ([49]) 
is probably somewhat underestimated. The value is also obtained using the average value 



of T(r] 77) of |T 



Let us now turn to the 0(p 6 ) results. Our main fit with m s /m = 24 gives 

Q = 19.6 [Eq. (20)] ; 21.5 [Eq. (|§] . (50) 

substantially lower than the values quoted in Eq. (0). On the other hand, Q is rather 
sensitive to the input value for m s /rh. In Fig. |3](a) we have shown how Q depends on the 
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o 




Figure 3: 
fit. 




m s /m m s /m 

(a) Q (b) m u /md as a function of m s /m but other input fixed as in the 



mam 



input values chosen for m s /m. Otherwise the inputs are as in our main fit. Fig. 0(b) shows 
similarly the dependence of m u /md on m s /rh. 

The position of the previous mentioned ellipsis can be described as the ratio of SU(3) 
breaking effects versus SU (2) isospin breaking using 



R = 



m, — m 



m d - m u 



(51) 



The standard value [4(| used p-uj mixing and combinations of the baryon masses to obtain 
R = 43.5 ± 2.2. Both of these inputs are now known to have rather large corrections 
of 0(m q ). p-Lo mixing is treated in [|6| and results for the baryon masses with earlier 



references can be found in ||7|. We obtain for this quantity 

R w 31 [Eq. (20)] ; 37 [Eq. ([!§] 

One of the main results from Table [I] is 

Tn 

— = 0.46 ±0.09. 



(52) 



(53) 



In the previous estimate the error is increased w.r.t. the one of the fit in order to include 
two other sources of errors, about 0.02 from the error on the electromagnetic correction 
due to m 2 K± — m? K0 and another 0.02 from the rather large change that happened with the 
scale variation at // = 1 GeV. The other new result we obtain is the strong mass difference 



Sm^ = (m^i — m^o) 



QCD 



0.32 ±0.20 MeV. 



(54) 
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Figure 4: The ratio m u /rrid versus the inputs for L r ±, L r & the rest of parameters correspond 
to those of our main fit. 

The change with the result quoted in [|l| has two sources, about 0.04 MeV from the inclusion 
of the mass-corrections in the Kaon electromagnetic mass difference and the remainder 
comes from the 0(p e ) effects. 

A possibly more serious variation with input is due to the assumptions on L\ and L r & . 
Our tests of large N c as described in Sect. |5.2| are mainly in quantities dominated by vector 
exchange. We have therefore checked what happens when we vary the input value of L\ 
and Lg over a fairly wide range. The resulting values for m u /rrid are plotted in Fig. |j. It 
should be kept in mind that not all of the points shown have good fits to all inputs and 
most of the points actually fall in the range of Eq. ( |53|) . It is a rather clear conclusion that 
the value m u /irid = is far away from the numbers we have obtained. Unfortunately this 
means that the solution to the strong CP problem needs to be found elsewhere. 

10 Summary 

In this paper we have explained how isospin breaking effects can be incorporated at next- 
to-next-to-leading order in CHPT. The expectation that isospin breaking would not sig- 
nificantly alter our previous results has been confirmed and the results presented here are 
a first step in understanding how isospin breaking effects can be added consistently. 

We have reanalyzed the main fit at 0(p 6 ) of the CHPT LEC's in the presence of isospin 
breaking. The determination of the L\ is only marginally changed compared to the isospin 
symmetric analysis. As by product we obtain control on the quantity m u /rrid at next-to- 
next-to-leading order in terms of meson masses. This constitutes one of our main new 



19 



results. 

Furthermore using the preliminary new experimental results on we can obtain one 
of the Zweig Rule suppressed LEC, namely L\. It follows approximately a linear relation 
together with U % . 

We also revised the value of the isospin breaking quantities Q and R discussing the 
electromagnetic effects. In view of the low value of Eq. (|50"D a reanalysis of the 7] —>■ 3n 
decay at the same order in the quark mass expansion as was done here seems necessary. 
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